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ABSTRACT 



Aims. We have carried out a survey of 60 Herbig Ae/Be stars in the 3 micron wavelength region in search for the rare spectral features at 3.43 
and 3.53 micron. These features have been attributed to the presence of large, hot, hydrogen-terminated nanodiamonds. Only two Herbig Ae/Be 
stars, HD 97048 and Elias 3-1 are known to display both these features. 

Methods. We have obtained medium-resolution spectra (R ~ 2500) with the ESO near-IR instrument ISAAC in the 3.15-3.65 micron range. 
\^ ' Results. In our sample, no new examples of sources with prominent nanodiamond features in their 3 micron spectra were discovered. Less 
than 4% of the Herbig targets show the prominent emission features at 3.43 and/or 3.53 yum. Both features are detected in our spectrum of 
HD 97048. We confirm the detection of the 3.53 yum feature and the non-detection of the 3.43 yum feature in MWC 297. Furthermore, we report 
tentative 3.53 yum detections in V921 Sco, HD 163296 and T CrA. The sources which display the nanodiamond features are not exceptional in 
the group of Herbig stars with respect to disk properties, stellar characteristics, or disk and stellar activity. Moreover, the nanodiamond sources 
Q ' are very different from each other in terms of these parameters. We do not find evidence for a recent supernova in the vicinity of any of the 
^ nanodiamond sources. 

C/2 We have analyzed the PAH 3.3 yum feature and the Pfund(5 hydrogen emission line, two other spectral features which occur in the 3 micron 
' wavelength range. We reinforce the conclusion of previous authors that flared-disk systems display significantly more PAH emission than self- 
^ shadowed-disk sources. The Pf 6 line detection rate is higher in self-shadowed-disk sources than in the flared-disk systems. 
? Conclusions. We discuss the possible origin and paucity of the (nano)diamond features in Herbig stars. Different creation mechanisms have 
r> been proposed in the literature, amongst others in-situ and supernova-induced formation. Our data set is inconclusive in proving or disproving 
either formation mechanism. 
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1. Introduction 

Herbig Ae/Be (HAEBE) stars are pre-main-sequence objects 
of intermediate (1.5-8 Mq) mass. These sources exhibit large 
infrared (IR) excesses, for the greater part due to thermal emis- 
sion of circumstellar dust. The geometry of this dust suiTound- 
ing t he late-B, A and F-type members of this clas s is disklike 
(e.g. [Mannings & Sargentlll997t iTesti et alJl2003t IPietu et alJ 
l2003t iFuente et al."2003' ' Natta et alJl2o'o4r Also for ea rly-B- 
type stars, evidence for disks has been found (e.g. Vink et al.l 
l2002h . However, there is growing evidence that the structure of 
the circumstellar matter in early-type B stars is fundamentally 
different from Herbig Ae and T Tauri stars (e.g. .Fuente et alJ 
l2003UBik & Thil2004 . 



* Based on observations collected at the European Southern 
Observatory, La Silla, Chile (program numbers 071.C-0653, 072. C- 
0648, 073.C-0245 and 074.C-0248). 

** Postdoctoral researcher of the Research Fund KULeuven. 



The chemistry and mineralogy of the dust component of 
circumstellar disks around HAEBE stars has been studied with 
unprecedented precision based on near-to-mid-IR spectra (2- 
200 jum) provided by the Infrared Space Observatory (ISO, 
^Cessleret al. 1996). The launch of this satellite in 1995 was 
the start of a brand-new era for the research of protoplane- 
tary disks. The spectra revealed a large variety in dust prop- 
erties and species, from small carbonaceous molecules like 
polycyclic aromatic hydrocarbon molecules (PAHs) to silicate 
dust. Moreover, some sources were shown to contain crys- 
talline grains, similar to those found i n come ts in our own so- 
lar system (Waelkens et al. 1996; Malfait et alJll998al Il999t 
Ivan den Ancker et al. 2000a b; Meeus et al '2001'). All avail- 
able ISO l-lSjim spectra of HAEBE stars have been studied 
as a whole bv lXcke & van den Anckeil ( l2004l hereafter AV04). 

A handful of HAEBE so urces (HD 97048 and Elias 3-1, 
IVan Kerckhoven et alJ 12002 , and references th erein) and the 
post-AGB binary HR 4049 (iGebaUe et alJll989l) display pecu- 
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liar spectral features at 3.43 and 3.53 /^m. Comparison with 
laboratory spectra has convincingly shown that the carriers of 
these features are hydrogen-terminated nanometer-sized dia- 
monds (hereafter nanodiamonds, Guillois et al. 1999). Near- 
IR observations of HD 97048 by iHabart et all ( l2004bl) have 
spatially resolved the emission region of the strong 3.43 and 
3.53 fim features on a scale of 0.18" x 0.18". The obser- 
vations prove beyong doubt that the emission region is con- 
fined to the inner 15 AU of the circumstellar disk. More tenta- 
tive detections of the features were re ported in the literature 
(e.g. in HP 142527 and HD 100546. IWaelkens et alj[l993 
iMalfait et al.l 1 1998b', respectively), but only one additional 
3.53 jum feature was discovered persuasively (in MWC 297, 
iTerada et alJ^Oli) . Summing up, only two HAEBE stars dis- 
play both the 3.43 and the 3.53 yum feature, and one HAEBE 
object shows the 3.53 jjm feature. Considering that the ISO 
sample contains 45 targets, only a minor fraction of the 
HAEBE stars appear to have such a spectacular 3 micron spec- 
trum (AV04). 

The currently operating successor to the ISO satellite, 
Spitzer, is not equipped with a spectrograph which can observe 
in the 3 micron region of the spectrum. Ground-based obser- 
vations are therefore the only alternative. With the present pa- 
per, we intend to enlarge the sample of HAEBE stars that are 
observed in the 3 micron region and possibly identify new tar- 
gets which display nanodiamond emission. Therefore, we have 
employed the near-IR ISAAC instrument (ESO-VLT) and fo- 
cused on the wavelength region around 3.3 (polycyclic aro- 
matic hydrocarbons, hereafter PAHs), 3.43 and 3.53 ^m. (nan- 
odiamonds). In Sect. 121 we present the sample and the data 
reduction procedure. The detected features are discussed in 
Sect. 13 

2. The data set 

2.1. The sample stars 

The sample of HAEBE stars compiled f or this survey has been 
based on Table 1 and 2 in the paper of 'The et al.1 lll994l) and 
on the Hst published by Malfait et al. ( 1998a). For 19 of the 60 
observed targets, including HD 97048, spectroscopic data in 
the 3 micron region are present in the ISO data archive (AV04). 
The overlap between the two samples guarantees the continuity 
between the two data sets. This survey enlarges the number of 
high-quality 3 micron spectra of Herbig stars available to us 
from 45 (ISO) to 82 (ISO-i-new observations). 

2.2. Observations and data reduction 

The observations in the 3 micron region of the spectrum were 
performed using the medium-resolution spectroscopic mode of 
the ESO-VLT instrument ISAAC' , mounted on UTl (ANTU). 
For each target two wavelength regions were recorded: the first 
is centered around 3.3 and the second around 3.5 /im. The 
wavelength coverage of each spectrum is about 0.25 yum. The 
width of the unresolved telluric lines in the spectra is a measure 

' http://www.eso.org/instruments/isaac/ 



for the resolution of the spectra. This value is approximately 
15 A, corresponding to a spectral resolution of 2300. The spec- 
tra are oversampled with a pixel-to-pixel wavelength step of 
2.5A. The width of the slit is I", needed to achieve the medium 
spectral resolution. The sky background subtraction is done by 
chopping/nodding with a throw of 20" in the direction of the 
sHt. 

The 3 micron region is heavily affected by atmospheric ab- 
sorption features. Standard stars (STDs) have been measured 
in order to correct for this effect. The straightforward data re- 
duction procedure consists of a flatfield correction, wavelength 
calibration and multiplication with a responsivity curve. The 
latter is deduced from the STD spectrum, scaled to the right 
flux values using an appropriate Kurucz model. It was noted 
that the wavelength calibration of the target and standard star 
is never exactly the same. The offset, a fraction of a pixel up to 
a few pixels large, leads to artefacts in the resulting spectrum. 
We have applied a relative shift to the extracted STD spectrum 
to correct for this effect, in an attempt to reduce the residuals 
of atmospheric features. 

During a few nights, no STD measurements were available 
to correct for the telluric line absorption. Some more nights, 
flatfield and/or wavelength calibration files were lacking^. In 
these cases, STD and calibration measurements from other 
nights were used. This increases the problem of wavelength 
calibration significantly, introducing second-order differences 
between the target and the STD spectrum. Furthermore, the 
assumption that the atmospheric conditions are approximately 
the same when observing target and standard star through the 
same airmass is not valid at all in these cases. We have tried to 
correct for the first order effect of this problem by applying an 
additional factor (~0.8-1.2) to the airmass ratio of target and 
standard star. This cosmetic surgery can of course not compen- 
sate for night-to-night differences in e.g. humidity. 

The 3.3 and 3.5 micron spectra were reduced indepen- 
dently. When comparing the flux levels in the overlap region, 
flux-level jumps occur for a few sources. This is probably due 
to slit losses: the slit width is only a fraction of the atmospheric 
seeing. Variations in seeing between the observations of stan- 
dard star and target make that a different fraction of the source's 
flux is captured. This effect translates to a multiplicative factor 
in the resulting spectrum. Other possible causes for the flux- 
level jumps include errors in background subtraction or dif- 
ferences in atmospheric conditions between the observations 
at 3.3 and 3.5 micron. We have corrected the jump between 
the two resulting spectra by applying a multiplicative factor to 
both spectra and scale them to the mean value in the overlap 
region. In most cases, this factor is of the order of 0.6-1.5. We 
note that, since the PAH 3.3 micron and nanodiamond 3.53 mi- 
cron feature are both included within one spectrum, their shape 
is insensitive to offsets in absolute flux level. The strength of 
these features is of course altered due to the application of the 
multiplicative factor. 



2 No STD in the 3.3 fim region on 04-05/09/2003 and 06- 
07/10/2004; no STD in both the 3.3 and 3.5 //m regions on 18- 
19/1 1/2003 and 17-18/07/2004; no calibration files on 04-05/1 1/2003 
and 20-21/10/2004 
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Fig. 1. The ISAAC spectra of the sample stars, rebiimed to a resolution of 750. When available, ISO-SWS {R = 150) or PHT-S 
(R ~ 90) spectra are overplotted in grey. The part of the ISAAC spectrum between 3.3 and 3.32 fim has been excluded. This part 
is unreliable due to the high atmospheric opacity in this wavelength range. 
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Fig.l. — (Continued) 
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Fig. 1. — (Continued) 
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Inherent to the reduction method, some residuals of tel- 
luric lines remain present in the final spectra. This is never- 
theless not dramatic for our survey. Because of the medium 
spectral resolution of our ISAAC spectra (~2300), the telluric- 
line residuals are very narrow compared to the typical width 
(FWHM~0.03//m) of the PAH and nanodiamond features. 
Furthermore a bright 3.53 yum source similar to HD 97048 
would be easily recognized, even without atmospheric correc- 
tions, due to the strength of the nanodiamond features. 

The reduced ISAAC 3 micron spectra of the sample sources 
are presented in Fig.Q] We have rebinned the spectra to a spec- 
tral resolution of 750 for the sake of clarity. The broad PAH and 
diamond features remain clearly visible. The sharp, unresolved 
Pfund 5 line on the other hand is spread out over 20 pixels and 
is therefore undetectable in most of the spectra shown in the 
figure. The smooth+rebin procedure makes that also the unre- 
solved telluric-line residuals are spread out. In regions which 
are rich in telluric absorption lines, this may cause problems: 
the lines are overlapping each other and form absorption bands 
in the rebinned spectrum. This induces "absorption features". It 
is clear that, of the two wavelength regions, the 3.3 ;um spectra 
suffer the most from telluric line absorption. The region be- 
tween 3.3 and 3.32 jum has been cut out. The measured flux 
in this region is close to zero in both science and standard-star 
observations, due to strong telluric absorption. Division by the 
STD spectrum induces a strong artefact in this wavelength in- 
terval. Note however that the PAH 3.3 /im feature can still be 
discerned, since the onset and peak position of its spectral sig- 
nature are located at shorter wavelengths. 

When available, the ISO spectrum is overplotted in Fig.[n 
The continuum flux levels of the ISAAC and ISO data are 
within a 30% error bar of the mean value of both for most tar- 
gets. The experimental error bar on the ratio of the ISAAC con- 
tinuum flux and the L-band magnitude — converted to the right 
units — is about 60%. The latter indicates the poor absolute flux 
calibration quality of the ISAAC spectra. 

2.3. Characterizing the spectral features 

Spectral features that appear in the 3 micron spectra are in this 
analysis characterized with a few parameters. The integrated 
line flux (in W m"^), full width at half maximum (FWHM, in 
yum), the peak flux (in Jy) and peak position (in fim) are de- 
termined. We have estimated the statistical error on these mea- 
surements based on the noise on the spectra. The error on the 
line fluxes due to uncertainties on the absolute flux level cal- 
ibration is estimated to be 60%. This error is not included in 
the figures and tables presented in this paper Note that the 
shape and the peak-over-continuum ratio of the detected fea- 
tures is not affected by the absolute calibration of the spectra. 
In Sect.|3lwe discuss the observed spectral features and include 
a comparison with the ISO measurements described in AV04. 

2.4. Classifying the spectral energy distributions 

We ha ve classified the sample targets into the two lMeeus et alJ 
(I2OOII) groups based on the shape of the spectral energy dis- 



tribution (SED). This classification roughly represents the cir- 
cumstellar disk's geometry: group I sources have a strong 
mid-IR excess which is believed to originate from a flared 
outer disk. Group II objects have a sma ller mid-IR excess 
and non-flaring outer disks jPuUemond l2002lb . The sam- 
ple targets in this paper are classified based on Fig. |2l (see 
also van Boekel et al. 2003; Dullemond et al. 2003; Acke et aU 
I2OO4 AV04). In this dia gram, the ratio of near-IR and mid- 
IR luminosity is plotted versus the non-color-corrected IRAS 
[12]-[60] color Group I sources and embedded objects end up 
in the lower right corner of the diagram, while group II sources 
have higher near-IR to mid-IR luminosity ratios. 

Sources with an Ay value exceeding 4 have been classi- 
fied as embedded sources. The link between SED shape and 
circumstellar geometry is less clear in these cases. In Fig. |2 
we have also included another possible classification criterium. 
Almost all group I and all embedded sources have a 60 micron 
excess exceeding 10 mag. In group II most objects have an ex- 
cess smaller than 10 mag. Both criteria — the diagram and the 
60 micron excess — hence lead to more or less the same clas- 
sification. In Table|5]we have indicated the classification based 
on the diagram. 

Two sample stars, T CrA and T Ori, cannot be classi- 
fied due to a lack of IR photometric measurements. However, 
both sources display photometric variability: their optical light 
curves are characterized by Algol-like minima, commonly re- 
ferred to as UX Orionis (UXOR) behaviour. These minima are 
attributed to extinction events which are caused by dust and 
gas clumps moving in and out the line of sight ( Wenzel 1968; 
'Grinin et alJ Il988l: IVoshchinnikovlll990l: lOrinin et al...l998.) . 
Dullemon d et alJ ( l2003l) have argued that the clumps spring up 
from the puff'ed-up inner rim, and that all UXORs have self- 
shadowed disks. In other words, all UXOR stars should be 
group II members. The authors indeed find this correlation in 
a large sample of Herbig stars. Based on this conclusion, we 
classify the UXOR sample targets T CrA and T Ori as group II 
members. 

In the present sample, 15 sources have been classified as 
group I members, 32 as group II sources and 6 as embedded 
objects. 7 sources remain unclassified due to a lack of photo- 
metric measurements and/or upper limits. 

3. Description of the detected features 

In Table |5l we summarize the detection of features in the 3 mi- 
cron spectra of the sample stars. A detection is defined as being 
a measurement of a spectral feature with a line flux exceeding 
3cr. Tentative detections are features with line fluxes between 
1 and 3cr. We have searched for the narrow Pfund 6 line at 
3.296 yum, the PAH 3.3 yum feature and the nanodiamond fea- 
tures at 3.43 and 3.53 yum. In five sources, also the Humphreys 
series of the hydrogen atom was detected. 

3.1. The Pfund 5 emission line 

The Pfund 5 line at 3.296 yum is located at the longer wave- 
length end of a region which is rich in telluric lines. The line 
is spectrally unresolved in all cases (i.e. FWHM of the order of 
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Fig. 2. The classification diagram of the sample targets. 
Sources with Av values above 4 are definied to be embedded 
sources (grey squares). Objects on the left side of the empirical 
separation line (dashed line) are group II members, on the right 
side group I. A classification solely based on the 60 micron 
excess gives almost equal results. Objects with a 60 micron ex- 
cess exceeding 10 mag (black diamonds) are almost all group I 
sources, while targets with an excess below 10 mag (grey di- 
amonds) belong mostly to group II. Some targets have been 
marked with their sequence number (see Table|5}; those are the 
(tentative) nanodiamond sources (see Sect. l3.4> . 



15A or equivalently 130 km s '). We have determined the char- 
acterizing parameters of the Pf S line at the full resolution of the 
ISAAC spectra (R=2300). The error bars on the line flux are 
consequently rather large. They are calculated from the noise 
in the 3.3 micron region of the reduced spectrum. The error in- 
cludes the photon noise, as well as the noise induced by telluric 
line residuals. The latter is by far more important. 

In the ISAAC sample of 60 sources, we detected the 
Pf 6 emission line in 16 targets or 27%. Four more objects ten- 
tatively display this feature, which adds up to a detection rate 
of 33% in this sample. 

The Pfund^ line has been (tentatively) detected in 
only three group I sources (HD 259431, HD 100546 and 
HD 142527). The detection rate in this group hence is 13-20%. 
In group II, 10-13 sources (tentatively) display the feature 
(31-41%), while half of the six embedded sample stars show 
the emission line (50%). In the group of unclassified sample 
stars, the detection rate is 14%. The Pfund^ line is more de- 
tected in group II and in the embedded sources, than in group I. 
This difference is not statistically significant on a 5% level. 

3.2. The Humphreys series 

Five of the sample stars (8%) show the Humphreys series of 
the Hi atom in emission: GU CMa, HD 76534, HP 130437 
V921 Sco and MWC 297 (for the latter, see also lTerada elai] 
l200lh . These objects are aU early-type stars (08-B2). The SED 
of GU CMa, HD 76534 and HD 130437 does not show a signif- 
icant near-IR excess. This indicates that these sources are prob- 



ably normal Be stars, rather than Herbig Be stars. The SEDs 
of V921 Sco and MWC 297 on the other hand are dominated 
by thermal emission of the circumstellar dust. The extinction 
in the line of sight towards these targets is strong. Ay - 5.1 
and 7.8 respectively. These massive stars are most likely sur- 
rounded by accretion disks. In Fig.|3] the 3.5 micron spectra of 
GU CMa, HD 76534, HD 130437, V921 Sco and MWC 297 
are presented. The Humphreys series is clearly visible. 

The presence of the hydrogen Humphreys series in the 
3 micron spectra correlates well with the presence of the hy- 
drogen Pfund 6 line: the Pf 6 line is detected in the spectra of 
GU CMa, HD 130437, V921 Sco and MWC 297. The 3.3 mi- 
cron spectrum of HD 76534 is rather noisy. We have not de- 
tected an emission line at 3.296 fim exceeding the Icr level in 
the latter case. 

3.3. The PAH 3.3 micron feature 

The 3.3 jum feature, attributed to the streching vibrations of the 
CH bond in PAH molecules, is detected in 7 Herbig sources 
(12%). Four additional sources tentatively show the feature. 
This suggests an upper limit for the detection rate of 18% in 
the ISAAC sample. The emission band has previously been 
detected in the infrared spectra of sample stars HD 97048, 
HD 100546 and HD 179218 (see e.g. AV04). We report new 
detections for HD 34282, V599 Ori and AS 310. Tentative de- 
tections are obtained for V346 Ori, V590 Mon, HD 101412 
and HBC 687. Also the ISAAC and ISO spectra of TY CrA 
contain the 3.3 micron feature. Nevertheless, this source is 
most likely not a HAEBE star, considering the absence of Ha 
emission and the fact that the IR emission emanates from the 

"TY CrA bar" which is n ot directly associated to the object 
I 1 1 

ISiebenmo rgen et al. 2000). For HD 34282, we have previously 
reported a non-detection of the PAH 3.3 micron feature in the 
ISO spectrum (AV04). The upper limit on the line flux de- 
duced from the ISO spectrum is 28.1 x 10"'^ Wm"^, which 
is consistent with the current detection in the ISAAC spec- 
trum (24 + 7 X 10""' Wm-^). The detections of the PAH fea- 
ture in the ISO spectra of HD 142666, HD 142527, V921 Sco 
and MWC 297 are not confirmed in the present spectra. For 
HD 142527 and MWC 297, strong Pf6 emission is observed, 
which is probably misinterpreted as PAH emission in the - 
lower spectral resolution- ISO spectra. Also in V921 Sco 
(=CD-42°1 1721), strong Pf6 emission is detected. The mea- 
sured line flux in the ISAAC spectrum is however not sufficient 
to explain the measured ISO 3.3 fim line flux. Remarkably, also 
in the ground-based TIMMI2^ slit spectrum of V921 Sco in 
the 10 micron wavelength region, no PAH features are detected 
(Fig. 1^. This may be surprizing, since the ISO spectrum of 
this source is dominated by strong PAH emission. The expla- 
nation is found in the extent of the emission region. Because of 
the large aperture size (14" x 20"), the ISO-SWS instrument 
has picked up PAH emission from the wider surroundings of 
the targets. Slit spectroscopy only selects a much narrower re- 
gion around the central source. For V921 Sco, this remark is 
in agreement with the result of AV04. Table 9 in that paper re- 
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Fig. 3. The ISAAC 3.5 micron spectra of GU CMa, HD 76534, HD 130437, V921 Sco and MWC 297. The Humphreys series of 
the Hi atom are indicated with vertical lines. The spectra in the top panels are clearly different from the spectra in the bottom pan- 
els. The first 3 targets are most likely normal Be stars, with the 3 micron continuum flux emanating from the stellar photosphere. 
The latter two are early-type accretion-disk systems. The near-lR continuum is due to hot circumstellar dust. 



veals that the luminosity of the observed PAH emission in this 
source roughly corresponds to what one may expect from a ten- 
uous halo. No emission features are detected in the 3.3 micron 
region of the ISAAC spectrum of HD 142666, possibly indi- 
cating that also in this case, the PAH emission is extended. 

In Fig.|5]we have compared the peak-over-continuum ratio 
of the PAH 3.3 micron feature measured in the ISO and ISAAC 
spectrum. The measurements/upper limits of HD 34282 (#2), 
HD 100546 (#34) and HD 179218 (#58) agree within the er- 
ror bai-s. For HD 142527 (#43) and MWC 297 (#50), sti'ong 
Pfund 6 emission led to a false PAH 3.3 jim detection in the ISO 
spectrum. In the previous paragraph we argued that the PAH 
emission in HD 142666 (#42) and V921 Sco (#47) is Hkely 
extended, and therefore not detected in the ISAAC spectrum. 
The ISO spectra of HD 97048 show significantly stronger PAH 
3.3, NANO 3.4 and NANO 3.5 ^lm features than the ISAAC 
spectrum. The relative ratios between the ISO and ISAAC line 
fluxes are 2.9+0.7, 2.4+0.2 and 1.65+0.02 respectively. The 
ISO aperture is many times larger than the slit width of ISAAC. 
The PAH and nanodiamond emission emanates from the cir- 
cumstellar disk's surface. The estimated extent is of the order of 
a few tens of an arc second ( van Roekel et al .'2004: Habart etTll 
l2004a.b.) . This extent, smeared out by the atmospheric seeing, 
could lead to a loss of PAH and NANO emission in the ISAAC 



spectrum, relative to the ISO spectrum. The underlying con- 
tinuum flux, which originates from the very inner parts of the 
disk, is less sensitive to this effect. This potentially explains the 
weakness of the features in the ISAAC spectrum. Alternatively, 
the PAH and NANO features might display temporal variabil- 
ity. This variability has not been noted in the five ISO 3 micron 
spectra which are available for this source, however These ISO 
observations span about one year and a half (AV04). 

The ISAAC detection rate of the PAH 3.3 feature in the 
group I sources is 36% (50% including tentative detections). 
The rate in group II is much lower (0-3% respectively). A 
PAH 3.3 yum feature was detected in the spectrum of the em- 
bedded source V599 Ori (17%). In the computation of these 
numbers, we have excluded GU CMa, HD 76534, HD 130437 
and TY CrA, which are most likely no Herbig stars. The high 
and very low detection rate in respectively group I and II con- 
firms the main conclusion of AV04. The PAH features emanate 
from the surface of a flared circumstellar disk. Only group I 
sources have a flar ed disk and are expected to display strong 
PAH features (Hab art et all20 04a. AV04). In Table[T]we com- 
pare the detection rates of the PAH 3.3 micron feature in the 
ISO sample to the rates in the present ISAAC sample. The 
ISO sample contains 16 group I, 22 group II and 7 embed- 
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Fig.4. The infrared spectrum of V921 Sco (CD-41°1 1721). 
The black curve is the ISO spectrum (R=150). The grey curve 
around 10 micron is the long-slit T1MM12 spectrum (cour- 
tesy of R. van Boekel). This spectrum is scaled to the IRAS 
12 micron flux. For clarity, a blow-up of the 3 micron wave- 
length region is shown in the inset. The grey curve in this panel 
is the ISAAC spectrum (R=750, not scaled); the black curve 
is the ISO spectrum at a resolution of 500. The slope of the 
ISAAC spectrum clearly deviates from the slope of the ISO 
spectrum at 3 micron, most likely due to slit losses. The ISO 
spectrum is dominated by strong PAH features. No PAH emis- 
sion is detected in the ground-based spectra. The ISAAC spec- 
trum shows the Pfundd line, the Humphreys series and pos- 
sibly a 3.5 jum feature. In the high-resolution ISO spectrum, 
the Pf5 line is visible on top of the PAH 3.3 yum feature. The 
TIMMI2 spectrum reveals that beneath the overwhelming PAH 
emission in the ISO spectrum, a shallow silicate absorption 
band resides. 
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Fig. 5. Comparison of the peak-over-continuum ratio of the 
PAH 3.3 jum feature, measured in the ISO and ISAAC spec- 
trum. Filled symbols indicate targets for which we have de- 
tected the feature in at least one of the infrared spectra. The 
numbers refer to the target numbers listed in Table |5] Upper 
limits are indicated by arrows. The dashed line represents ratio 
ISAAC = ratio ISO. 



Table 1. The detection rates of the PAH 3.3 //m feature in the 
ISO and ISAAC samples. T stands for the total of classified 
sample sources, I and II refer to the Meeus groups and E rep- 
resents the embedded sources. The last row gives the detection 
rates in the combined ISO+ISAAC sample. When two percent- 
ages are given, the second refers to the detection rate includ- 
ing tentative detections. The detection rate in group I is sig- 
nificantly higher than in group II (confidence level > 99.5%). 







T 




I 


II 


E 


ISO 


22- 


-29% 


44- 


-56% 


14% 


0-14% 


ISAAC 


12- 


-18% 


36- 


-50% 


0-3% 


17% 


total 


14- 


-21% 


32- 


-48% 


5-7% 


9% 



ded sources . The combined ISO-i-ISAAC sample consists of 
25 group 1, 42 group II, 1 1 embedded and 4 unclassified targets, 
in total 82 sources. A number of targets are present in both the 
ISO and ISAAC sample. The sources with a PAH detection in 
the ISO spectrum but not in the ISAAC spectrum, or vice versa, 
were all discussed in the previous paragraphs. We have consis- 
tently included these results in the computation of the detection 
rates. Up to half of the sources classified in group I display the 
PAH 3.3 jim feature, while it is detected in only a small number 
of group II and embedded objects. 

According to AV04, the PAH 3.3 micron feature is the least 
detected PAH feature in the 3-11 micron region in the ISO sam- 
ple. It is present in only 44% of the cases where the 6.2 micron 
feature is detected. The latter is the most detected PAH fea- 
ture in the ISO study. Furthermore, the PAH 6.2 feature is al- 
ways detected when the PAH 3.3 jum feature is present^. Using 
Bayes's theorem^, one can estimate the hypothetical detection 
probability of the 6.2 micron feature in the total sample. We 
assume that the detection probability of the 3.3 micron fea- 
ture is equal to the relative detection rate in the total sample. 
The conditional probabilities -the relative detection rate of one 
feature, given the presence of the other- are taken from the 
ISO study (AV04). For group I, the hypothetical PAH 6.2 mi- 
cron detection probability is 73%- 100%, while it is 1 1%-16% 
for group II and 20% for the embedded sources respectively. 
Hypothetically, three quarters to all flared disk sources are PAH 
emitters. Most of the PAH detections in group II and in the em- 
bedded sources, as well as a fraction of the group I detections, 
are likely due to extended emission from a circumsystemic neb- 
ula. 

3.4. The nanodiamond features 

It is clear that our ISAAC sample does not contain any targets 
with a high 3.53 yum feature/continuum ratio except HD 97048. 

* We have classified T CrA as group II member and V921 Sco as 
embedded source in stead of in group I (as was done by AV04). 

^ Except for HD 142527, a PAH 3.3 fim detection which was how- 
ever not confirmed by the ISAAC spectrum. 

P{B) = ^"^Lllio^''^' with A and B stochastic events. In this case A 

^ ^ P{A\B) 

and B represent the detection of the 3.3 and the 6.2 jim feature in the 
spectrum, respectively. 
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This in itself is an important result. The detection of a strong 
3.53 yum feature in the near-lR spectra of Herbig stars is thus 
Umited to Elias 3-1, HD 97048 and MWC 297. This leads to 
a detection rate of 3/82 - 4% for the entire ISO+ISAAC sam- 
ple. The 3.43 jim feature is even more rare, since it has been 
detected in Elias 3-1 and HD 97048 only (2% of the sample). 

We confirm the detection of the 3.53 micron feature in 
MWC 297 llTerada etalJbOOll AV04). The shape of this com- 
plex band resembles that of the 3.53 jjm feature in HD 97048, 
in the sense that the red wing is steeper than the blue wing 
and the peak position is around 3.529+0.002 yum. It is there- 
fore likely that these features have a common carrier, i.e. 
(nano)diamonds. The line flux ratio of the 3.43 and 3.53 ;um 
features is 0.82 (ISO) for EUas 3-1 and 0.41 (ISAAC; ISO: 
0.59) for HD 97048. The upper Umit to this ratio is 0.39 (ISO; 
ISAAC: 2.0) for MWC 297. This indicates that the NANO 
3.43 yum is, if present, at least weaker than in the two other 
(nano)diamond sources. 

A few more sample targets display a feature around 3.5 //m. 
In the Appendix, we argue that most of these features are arte- 
facts. Only for V921 Sco, HD 163296 and T CrA, a potentially 
real feature is detected. These detections are by far not sure, 
however. In the following, we refer to these sources as the can- 
didate nanodiamond emitters. 

4. What makes the nanodiamond emitting sources 
so special? 

It is clear that very few targets display nanodiamond features in 
their 3 micron spectra. The main question which is addressed in 
this paper concerns the origin of these features. In this section, 
we will show that the 3.43 and 3.53 jjm emitting sources are 
no outliers in the group of Herbig stars as a whole (1), but also 
that these targets have very diff'erent stellar and circumstellar 
parameters compared to each other (2). 

4.1. Nanodiamonds and polycyclic aromatic 
hydrocarbons 

PAHs and nanodiamonds are both carbonaceous grains. The 
3.3, 3.43 and 3.53 features have been attributed to the C-H 
bonds in these structures. The emission of both species has 
been shown to be spatiall y extended and linked to the surface 
of the circumstellar disk ( v an Boekel et al. 2004: Habar t et al.1 



O 100.0 



(D 

C 10.0 



l2004a.b.. AV04). Fig.|6ldisplays the PAH 3.3 jum Hne flux ver- 
sus the nanodiamond 3.53 yum line flux. The line fluxes deter- 
mined for Elias 3-1 are included as well. The latter are based 
on the ISO spectrum of this source, rebinned to a resolution 
of 500. For the sake of clarity, we left out the upper limits 
when both features were undetected. Although a correlation 
between PAHs and nanodiamonds may be expected based on 
their common chemical root, the PAH 3.3 fim feature is not de- 
tected in MWC 297 and the three tentative nanodiamond emit- 
ters V921 Sco, HD 163296 and T CrA. Furthermore, the PAH 
3.3 jum feature is obviously much more common than the nan- 
odiamond features in the group of HAEBEs. There does not 
seem to be a connection between them. However, it should be 
noted that the line flux ratio of PAH 3.3 and NANO 3.53 jum 
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Fig. 6. The line flux ratio of the PAH 3.3 jum and the nanodi- 
amond 3.53 yum feature versus the NANO 3.53 fim line flux. 
The dashed line indicates where both line fluxes are equal. 
The black square indicates the ISO measurements of Elias 3-1. 
The black diamonds represent HD 97048 and MWC 297. The 
grey diamonds are the nanodiamond candidates V921 Sco, 
HD 163296 and T CrA. The black dots represent the targets 
with a detected PAH 3.3, and undetected NANO 3.53 feature. 
Only in HD 97048 and Elias 3-1, both features are detected. 
The line flux ratio of the PAH 3.3 and NANO 3.53 features is 
smaller than one in these targets, while the upper limits on this 
ratio are larger than unity in the other (candidate) nanodiamond 
emitters. 



in HD 97048 and Elias 3-1 -the only two sources where both 
features are detected- is significantly smaller than 1 (0.13+0.03 
and 0.5+0.1, respectively). The upper limit for this ratio in the 
other four (1+3 candidate) nanodiamond targets is larger than 
unity. Hence it may be that the PAH 3.3 /urn feature is present, 
but too weak to be detected in the latter sources. If it is always 
true that the line flux ratio of PAH 3.3 iJ.m and NANO 3.53 jum 
is smaller than one, the upper limits on the NANO 3.5 line flux 
in the other stars with detected PAH emission indicate that no 
nanodiamond feature can be present; the lower limits on the 
line flux ratios are all above one in these targets. In conclusion, 
the nanodiamond feature at 3.53 fim appears to be stronger than 
the PAH 3.3 fim feature. In the other sense, the presence of a 
PAH 3.3 fim feature does not imply the presence/absence of a 
nanodiamond feature. 

4.2. Appearance of the disk 

As mentioned before, the (nano)diamond features are linked to 
the circumstellar disk. Here we investigate whether the disks of 
MWC 297, Elias 3-1 and HD 97048 are in any way different 
from other HAEBE disks. We have summarized some over- 
all disk parameters of the (candidate) nanodiamond sources in 
Table|2 As a reference, we list the median value of each param- 
eter in each group. The median values of Lgxc/i^* and the spec- 
tral index (a proxy for dust particle size) are not significantly 
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Table 2. Some disk parameters of the (candidate) nanodia- 
mond sources. The Meeus group indicates the overall geome- 
try of the disk. The ratio of the total IR excess to the stellar 
luminosity, Lexc/i^*, indicates the fraction of stellar flux that is 
captured and reprocessed by the disk. For passive disks, this 
number represents the fraction of the sky covered by the disk, 
as seen from the central star. The ratio for T CrA is likely 
overestimated due to strong nearby IR emission, not directly 
linked to the system. The parameter n is the spectral index in 
the 450 jum-1.3 mm range, defined by AF^ oc A^". 



# Target 


Group 


^exc 1 l^-k 


n 


Elias 3-1 


1 


0.5-1.3" 


3.23±0.03 


33 HD 97048 


1 


0.40 


>2.9 


50 MWC 297 


E 


0.06 


3.1±0.2 


47 V921 Sco 


E 


0.11 


4.1±0.3 


49 HD 163296 


11 


0.34 


2.9±0.1 


57 T CrA 


11? 


3.8 


7 


Median 


1 


0.50 


3.9 


Median 


11 


0.44 


3.0 


Median 


E 


0.63 


4.1 



° For Elias 3-1 we used a stellar luminosity of 20-50 Lq 
Table|3j. 



(see 



different in sources with nanodiamond features as opposed to 
the sources without. 



4.3. Stellar characteristics 

MWC 297, EUas 3-1, HD 97048 and the three candidate nan- 
odiamond stars have very different stellar parameters. We have 
summarized their spectral type, estimated age and stellar lu- 
minosity in Table 13 The median parameters for the total sam- 
ple are again given as a reference. Nanodiamond features are 
found in sources with a wide range of spectral types (B1.5- 
FO). Although the nanodiamond sources might be a bit older 
(> 1 Myr) than typical for Herbig stars, there exists no strong 
correlation with age; there are many examples of Herbig stars 
with similar ages in our sample which do not show nanodia- 
mond features of comparable strength. 

The absolute strength of the nanodiamond features (i.e., the 
nanodiamond luminosity, included in Table appears to be 
linked to the spectral type (and hence Tgff and luminosity) of 
the central star. This weak correlation appears to be the only 
one which can be derived. The nanodiamonds in the circum- 
stellar disk must be heated to a certain temperature to be able 
to radiate in the near-IR. Hotter stars heat larger volumes of 
circumstellar dust. The extent of the nanodiamond emission re- 
gion is therefore expected to be larger around earlier-type stars. 
This is in agreement with the observed correlation. 

One of the current hypotheses on the presence of the nan- 
odiamonds in our own solar system suggests that they were 
formed in a supernova outflow, and injected in to the system 
(idea originally put forward by IClavtonllTgsi see Sect. Bl. 
Under this hypothesis, the presolar nebula is enriched in mat- 
ter from the supernova outflow. The parent star of a disk which 
contains supernova-injected diamonds must therefore be pol- 



Table 3. Some stellar parameters of the (candidate) nan- 
odiamond targets, and the observed nanodiamond 3.53 fim 
luminosity Lag. The median value of each stelllar parame- 
ter in the total sample is given as a reference. Values taken 
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# 


Target 


Sp. Type loj 


i Age/yr loj 








Elias 3-1 


A6:e 6 


1.3-6.9 1.3-1.7 


1 X lO^' 


33 


HD 97048 


B9.5Ve+sh 


>6.3 


1.9 


5 X 10-2 


50 MWC 297 


B1.5Ve 


5-6 


4.2 


1 X 10-' 


47 


V921 Sco 


BOlVep 


5-6 


4.9 


8 X 10-- 


49 


HD 163296 


A3Ve 


6.6 


1.6 


1 X 10-3 


57 


TCrA 


FOe 


7 


0.8 


3 X 10-* 




Median 


A2 


6.0 


1.7 





luted by supernova material as well. Supernova outflows har- 
bor r- and p-process elements and specific isotopes of lighter 
elements. Enrichment of the stellar photosphere by such an out- 
flow is expected to alter the photospheric abundance pattern of 
the star significantly. An abundance analysis of the nanodia- 
mond sources opposed to stars with comparable stellar param- 
eters might clarify whether or not supernovae are indispens- 
able for the formation of (nano)diamonds in their circumstellar 
disks. For HD 97048, two high-quality optical spectra are avail- 
able to us: a FEROS 3900-9100A spectrum with a signal-to- 
noise ratio (S/N) of 115, and a UVES 4800-6800A spectrum 
with a S/N of 460. The latter spectrum is compiled from 12 
UVES spectra, obtained in a time spa n of a fe w days. Th e spec- 
tra have been previously published i n lAcke & Waelken1( E004l) 
and lAcke & van den Anc ker ( 2006), respectively. Based on the 
Vienna Atomic Line Database (VALD, Kupka et aL1999.) . we 
have looked for the strongest absorption lines of elements with 
atomic number larger than 30 (Zn). At none of the VALD wave- 
lengths, an absorption line is detected. We have deduced upper 
limits on the equivalent width of the strongest lines and hence 
on the abundances of four heavy elements. The abundances of 
Y, Zr and Ba in HD 97048 are less than a few times the so- 
lar abundance. Strontium is even depleted by at least a factor 
10 compared to our sun. We conclude that there is no observa- 
tional evidence for a photosphere enriched with heavy elements 
in the nanodiamond source HD 97048, and hence no evidence 
for the occurance of a recent supernova event in the target's 
proximity. 

4.4. Disk and stellar activity 

In Fig. the line flux ratio of the Pf 6 line and the nanodi- 
amond 3.53 fj.m feature is plotted versus the NANO 3.53 jum 
line flux. Group I members HD 97048 and Elias 3-1 have no 
detected Pf 6 line. The massive young stellar object MWC 297 
does display both features. We have indicated the measure- 
ments for the tentative nanodiamond emitters (V921 Sco, 
HD 163296 and T CrA) with grey plotting symbols. Also for 
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Fig. 7. The line flux ratio of Pf 5 and NANO 3.53 ;um versus 
the NANO 3.53 fim line flux. Legend see Fig.|6l 



these three targets, there is no uniformity. The early-type em- 
bedded source V921 Sco and group II member HD 163296 dis- 
play the Pf S line, while the UXOR object T CrA does not. 

X-ray emission in stars with convective envelopes, like the 
lower-mass T Tauri stars, is often explained as being due to the 
presence of a magnetic field. The matter is less clear for the 
earlier- type Herbig stars. Although these stars are believed to 
have radiative envelopes and thus cannot drive a solar-like mag- 
netic dynamo, the X-ray luminosity of Herbig stars is neverthe- 
less explained in the ligh t of magnetic effects. Such a magnetic 
field can be fossile (e.g. ISkinner et al .' 2004: H amaguchi et al.l 
or produced differently than in our Sun (e.g. the differ- 
ential rotation dynamo model. Tout & Pringle 1995). Hard and 
soft X-ray spectra are observed in the group of HAEBE stars. 
Generally, the first kind is attributed to hot plasma in a stellar 
corona or the presence of a late-type companion, the second 
to accretion events onto the photosphere. In O-type stars, ener- 
getic wind shocks can produce X-ray emission. It appears that 
Her big winds are not powerful enough to produc e the same ef- 
fect dSkinner et al.l2004tlHamaguchi et alJ2005l) . 

Whatever mechanism causes the X-ray luminosity in 
Herbig stars, it does not seem to influence the presence/absence 
of the nanodiamond features. In Fig. |8] we have plotted 
the X-ray luminosity versus the NANO 3.5 l ine lurn inosity. 
The X-ray luminosities are taken from Skin ner et alJ §004), 
iHamaguchi et alJ (l2005l) and IPamiani et al.l J2006I) . It appears 
that, in this small sample of observed Herbig stars, the (candi- 
date) nanodiamond emitters have typical X-ray luminosities. 
When rescaling the NANO 3.5 line flux to the appropriate 
distance to the source, MWC 297 appears to be the most lu- 
minous nanodiamond target (La g - 0.10 L©). It is a rather 
strong X-ray source, although the X-rays likely emanate from 
nearby lower-mass stars (we follow the suggestion that 5% or 
less of the total X-ray l uminosity around MWC 297 enianate s 
from the B 1 .5 star itself. lVink et a l."2005':'Dami ani et al.l20 06). 
Elias 3-1 has a comparable X-ray luminosity, but a 3.53 //m lu- 
minosity which is ten times smaller. HD 97048, the third con- 
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Fig. 8. The X-ray luminosity Lx versus the NANO 3.53 jim 
line luminosity L3 5. The latter quantity is the line flux rescaled 
by a fac tor An cfi, w it h d t h e distance to the source as 
hsted in IS kinner et al.' ('2004|), IHamaguchi et alJ ll2005l) or 
IPamianiet al..c2006 ). We apply the same distance factor to the 
3.53 /im line flux as was applied for the cited Lx, to avoid sys- 
tematic errors. The plotting symbols are the same as in Fig.|6l 
The (candidate) nanodiamond emitting sources have X-ray lu- 
minosities which are not uncommon in the group of sample 
stars observed in the X-ray wavelength range. 



firmed nanodiamond emitter, has an intermediate nanodiamond 
luminosity (L3 5 = 0.05 L©), but an X-ray luminosity 15 times 
more modest than the previous two targets. Recently, Elias 3-1 
has been identified as a binary (with a separation of 9 AU at 
a distance of 160 pc. 'S mith et alJ l2005'). The companion is a 
lower-mass star, possibly with a convective envelope. The lat- 
ter may explain the rather high X-ray luminosity attributed to 
Elias 3-1. 

The analysis presented in the previous sections suggests 
that the few sources which display nanodiamond features 
are not exceptional compared to the bulk of Herbig stars. 
Furthermore, these targets do not appear to have obvious stellar 
or circumstellar characteristics in common. 



5. Conclusions and Discussion 

We summarize the two main conclusions of the present paper. 

- Our survey of 3 micron spectra of HAEBE stars has 
not revealed a new spectacular nanodiamond source like 
HD 97048 or Elias 3-1. This negative result increases 
the rarity of these near-IR bands in the group of Herbig 
stars. Within the current sample, only 4% of the targets 
are surrounded by detected nanodiamonds. We find a large 
range in nanodiamond luminosity, which spans 1 -2 orders 
of magnitude. No clear correlation is found between the 
presence of the nanodiamond features and the disk and 
stellar parameters, although a weak link between nanodi- 
amond luminosity and stellar luminosity may be present. 
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Furthermore, there does not appear to be a specific disk/star 
property which the nanodiamond sources have in common. 
We have shown that our findings are inconclusive with re- 
spect to the influence of the surroundings of the targets. We 
do not find any evidence for a correlation of nanodiamond 
emitters with the locations of known recent supernova rem- 
nants, nor do we find evidence for an enrichment in heavy 
elements of the photosphere of the system with the most 
prominent nano-diamond features, HD 97048. 
- Most to all Herbig stars with a flared disk are PAH emit- 
te rs. This conclusion once again confirms the suggestion 
ofl Meeuse'tal]ll200ll) that PAH molecules can only be ex- 
cited in systems where they can capture direct stellar UV 
photons. In the case of self-shadowed disks, the puffed-up 
inner rim casts a shadow over the outer disk, preventing 
the PAH molecules from being heated by the c entral star. 
Our find ings reinforce the results of AV04 and Habar t et"al] 
The latter studies were based on ISO data. Due 
to the large aperture sizes of the ISO instruments, it is diffi- 
cult to differentiate between PAH emission emanating from 
the disk surface, and contamination from the wider sur- 
roundings of the target. The present ISAAC data set, which 
clearly confirms the correlation, is considerably more re- 
liable in this sense: both the narrow slit and the chop-nod 
technique, which allows for an efficient subtraction of con- 
stant background emission, reduce the possibility of confu- 
sion with extended PAH emission drastically. 



Laboratory experiments have been instrumental in the 
identification of the carriers of the 3 micron features 



(e.g. Guillois et al. 


"19991 ISheu et al.1 120021: iTones et alJl2004 


IMutschke et al. 


.2Q0^. Currently, the astronomical spectra 



can be reproduced from the a bsorption spectr a of films of 
hydrogen-terminated diamonds ( Sheu et al. '2002V The spectral 
features in the 3 micron range are due to CH-bond stretches (as 
is the PAH 3.3 /zm feature). The prominent peaks around 3.43 
and 3.53 yum observed in HAEBE spectra can only be repro- 
duced when the diamonds are sufficiently large. Typical sizes 
for the diamonds in HAEBE stars (~50-100 nm) are 2 orders of 
magnitude larger than the nanodiamonds found in our solar sys- 
tem (e.g. Jones et al. 2004 ). The spectral features hence indi- 
cate the presence of large, hot hydrogen-terminated diamonds. 

The question remains whether there are indeed no di- 
amonds in the majority of disks. It may be possible that 
nanodiamonds of small sizes are present, but remain unde- 
tected, because th eir spectral signature is less pronounced (e.g. 
IShenetallEooi) . Also diamonds which are not hydrogen- 
terminated are not detectable in the 3 micron region. Diamonds 
that are stripped from hydrogen atoms are hence overlooked by 
our survey. A third possibility why (nano)diamonds can remain 
undetected is their exact location — and hence temperature — 
in the disk. Cool hydrogen-terminated diamonds will not radi- 
ate significantly in the near-IR. The absence of the 3 micron 
features in the spectrum of a HAEBE system hence does not 
necessarily mean that no (nano)diamonds are present. 



In-situ formation or extraneous origin? 

Diamon d was the first pre solar mineral to be identified in me- 
teorites jLewis et al.l 19871) and it has the highest relative abun- 
dance among carbonaceous grains. The diamond has been clas- 
sified as presolar based on the presence of a significant over- 
abundance of very heavy (Xe-H) and very light (Xe-L) Xe 
isotopes (together Xe-HL). The overabundance of this noble- 
gas carrier is expected to be produced by r- and p-process nu- 
cleosynthesis in supernovae. The link between the presence of 
nanodiamonds and the xenon overabundance led to the hypoth- 
esis that the diamonds are f ormed in, and in jected into the sys- 
tem by a supernova outflow llClavtonll989l) . Isotope anomalies 
on other heavy elements such as Te and Pd also point to super- 
nova nucleosynthesis (Richter et al. 1998: Maas et al, 20011) . 

Recent analysis of ^"Fe isotopes in chondrites has clearly 
demonstrated that a nearby type II supernova has contributed 
material to the natal cloud from which our own solar s ystem 
formed dMostefaoui e t al.1 1200.4 iTachibana et aLl l2006l) . This 
process may have been important for planet formation, as the 
the radioactive decay of ^"Fe into ^**Ni was an important heat 
source for the early planetary melting and dififerentiation and 
keeping asteroids thermally active for much longer than would 
be possible from the decay of ^^Al alone. 

Infrared spectroscopy has shown that solar-system dia- 
monds contain N and O, probably i n chemic al functional 
groups on their surfaces jLewis et al.li 1989,: Mutschke etalJ 
ll 99.4 lAndersen et al.lll 998t I.Tones et al Jl2004h . However, the 
similarity of C- and N-isotopes of these diamonds and the solar 
system as a whole supports the idea that not all solar-system di- 
amonds originate from supernovae and that the supernova con- 
tribution to the diamonds in our own solar system is probably 
not very large. It thus remains unclear how the majority of the 

solar-system diamonds has formed. 

Al though ad vocated by several recent studies ("Ko uchi et alJ 

l2005t Ib inette et al 2005.) . the presence of nanodiamonds in 

the interstellar medium remains controversial. The absence of 
observable 3.43 and 3.53 features in the ISM suggests 
that hydrogenated hydrocarbons cannot b e more abundant tha n 
0.1 ppm in the interstella r medium dTielens et al.l Il987h . 
IVan Kerckhoven et al.1 (|2002') showed that in HD 97048 and 
Elias 3-1,1 nanodiamond part per billion relative to hydrogen 
is required to reproduce the 3 micron spectra observed by ISO, 
making an ISM origin of nanodiamonds in principle a viable 
option. This hypothesis has problems in explaining the paucity 
of detections of the 3.43 and 3.53 fim feature in Herbig stars 
within our sample, however. 

We found no evidence for the presence of a supernova rem- 
nant (SNR) in the proximity of any of our sample sourc es. We 
made use of the catalog of galactic SNRs published bv lOreenl 
(2004). Only R Mon and V590 Mon are potentially located 
close to a SNR; the Monoceros Nebula (G205. 5+0.5). These 
sources display no detected nanodiamond features, however. 
We note that the absence of a SNR in the vicinity of our sample 
stars in the above-mentioned catalog does not necessarily im- 
ply that no supernova has actually occured. Most of the sample 
targets are located in massive-star-forming regions. It may very 
well be that the most massive members in such a region have 
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recently gone off as a supernova. In a relatively crowded star 
forming region, the ejecta of such an event likely affect a num- 
ber of surrounding young stellar objects or molecular clouds. 
In this scenario, the presence of diamonds in the circumstellar 
disk is expected to be common to all disk systems in the vicin- 
ity of the supernova. Unfortunately, we have no spectra of disk 
sources close to Elias 3-1, HD 97048 or MWC 297 to check 
this hypothesis. 

The parent star of a disk which contains supernova-injected 
diamonds must be be polluted by supernova material as well, 
under the supernova hypothesis. Enrichment of the stellar pho- 
tosphere by a supernova outflow is expected to alter the pho- 
tospheric abundance pattern of the star significantly. An abun- 
dance analysis of the strongest diamond source HD 97048 has 
produced upper limits on the abundances of a few heavy el- 
ements. The Sr abundance is at least one order of magnitude 
less than the solar abundance of this element. No evidence sup- 
porting the supernova hypothesis is found. 

Except for HD 97048 and Elias 3-1, both nanodiamond 
features h ave also been observed in the post-AGB binary 
HR 4049 jGebaUe et al.|[l?89i) . The post-AGB phase is very 
short (~ 10,000 yr) compared to the lifetime of an intermediate- 
mass pre-main-sequence star. Due to the short timescales, it 
is unlikely that a nearby supernova is the cause of the pres- 
ence of nanodiamonds in this system. However, the oxygen 
isotopes in the circumbinary disk of HR 4049 display pecu- 
liar behaviour: the relative '^O and '^O abundances are two 
orders of magnitude larger than the surface abundances in 
other evolved stars ('^O/'J O = 8.3+2.3 and "'O/^'^O = 6.9+0.9, 
ICami & Yamamurall200ll) . These exceptional values are most 
likely related to the binary nature of HR 4049, and possibly to 
a nova outburst. The presence of nanodiamonds in the circum- 
stellar environment may be related to such an event as well. 

The extraneous or igin of the circumstellar diamonds is dis- 
puted, however. iDai et all ll2002) have stated that, if the super- 
nova hypothesis is true, nanodiamonds should be abundant in 
solar-system comets as well. These objects are believed to have 
formed further out in the early solar system and are likely more 
pristine than meteorites. The authors have investigated fragile, 
carbon-rich IPD particles which enter the Earth's atmosphere 
with speeds in the range of cometary bodies, which suggests 
that the grains originate from these objects. Nanodiamonds are 
absent or strongly depleted in such IPD grains, which indicates 
that diamonds are not uniformly abundant in the solar system. 
This may support the hypothesis that the detected meteoritic 
nanodiamonds have formed in situ, and are not of presolar ori- 
giiij 

iGoto etal] ( l2000l) has suggested that the X-ray hardness of 
the radiation field of the central star may play a decisive role in 
the formation of diamond. However, our analysis has shown no 
obvious link between the source's X-ray luminosity and pres- 
ence/absence of nanodiamonds. Furthermore, the strength of 
the few detected diamond features appears to be independent 
from the X-ray strength of the central star as well. 

The extraneous origin of nanodiamonds is an attractive hy- 
pothesis, because it naturally explains the paucity of nanodi- 
amond sources and the seemingly very different properties of 
these targets. No evidence was found, however, to substanti- 



ate the claim of an external triggering source in the vicinity of 
the targets. In-situ formation of nanodiamonds is a viable al- 
ternative. However, it remains unclear why only very few disks 
can/could provide the specific conditions needed to produce di- 
amond. The disk and stellar parameters of the diamond sources 
are not uncommon, nor do the targets form a consistent group 
within the Herbig stars. Selection effects in terms of the size, 
temperature or hydrogenation of the diamonds may nonethe- 
less be present. Such effects could prevent the detection of fea- 
tures in the 3 micron region in spectra of the majority of cir- 
cumstellar disks. A larger sample of systems displaying dia- 
mond emission is needed to distinguish between in-situ forma- 
tion and an extraneous origin of the (nano)diamonds in Herbig 
Ae/Be stars and in our own solar system. 

Appendix 

Apart from HD 97048 and MWC 297, a few other targets dis- 
play a feature in the 3.5 micron region. Due to their low fea- 
ture/continuum ratio, however, it is difficult to discern them 
from data reduction artefacts. In this Appendix we will try to 
clarify which features are astrophysical, and which are not. 

In our opinion, the main source of data reduction artefacts 
is the correction procedure for atmospheric absorption lines. 
Assume that 5^ is the flux emitted by the standard star (STD), 
and T,\ the atmospheric optical depth at a certain airmass and at 
the time of the STD observation. Define the scientific target's 
flux as T;^ and the corresponding atmospheric optical depth 
as (Ta- We assume that the airmasses at the time of both ob- 
servations are the same. Note however that even in this case 
ta and CT) are not necessarily the same, since the atmospheric 
conditions (e.g. air density, humidity) may differ between the 
two observations. The STD's flux can be approximated by a 
photospheric model. If F^. is the measured STD flux, then 
e'"' = S A,modei/FA Can be estimated. This atmospheric correc- 
tion factor is then applied to the target measurement, which 
leads to an estimate of the target's flux, - {T,i e"""-') e'"' - 
T,i e^'^i^'^i* a! T,i- Slight differences between and cr^ intro- 
duce artefacts in the reduced spectrum of the science target. 
As a first order approach, we assume that the main difference 
in telluric extinction occurs due to a geometrical and/or atmo- 
spheric density effect. Such an effect is linear in t^. Hence 
g(T.i-tri) proportional to e'^'^\ with C a constant. The bottom 
two spectra in Fig.|9lshow the influence of over- and undercor- 
rection on a flat continuum. The artefacts that are relevant in 
this discussion are the "emission features" at 3.505 /im in the 
case of undercorrection (oc e"'"'), and 3.54 //m in the case of 
overcorrection (oc e^*). 

We will try to separate the real spectral features from the 
artefacts in the 14 candidate nanodiamond emitters listed in 
Tabled Two of them, HD 97048 and MWC 297, are known to 
display the 3.53 feature and will serve as a reference. The 
vast majority of the other 12 detected features is most likely due 
to the data reduction procedure. In Fig.|9l the normalized spec- 
tra of the candidate nanodiamond emitters have been displayed. 
We have indicated which spectra appear to be over- and under- 
corrected. Most of the features follow the spectral shape of e^' 
or e-^K Only for V921 Sco, HD 163296 and T CrA, die de- 
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tected feature's appearance deviates from the artefacts' shape. 
Even in these cases however, it remains difficult to exclude an 
artefactual origin. An indication may come from the upper limit 
to the 3.53 jum line fluxes measured in the ISO spectrum. The 
upper limits for HD 163296 and T CrA are in agreement with 
the measured ISAAC line fluxes (<2.5 and <1.4 vs. 2.3 and 
0.65 X 10"'^ Wm-2). For V921 Sco, the observed Une flux 
exceeds the ISO upper limit by a factor three (<3.7 vs. 10.3 
X 10"'^ Wm"^), which makes the NANO 3.53 yum detection 
in this source rather suspicious. In the paper, V921 Sco, HD 
163296 and T CrA are nevertheless treated as the remaining 
three candidate nanodiamond emitters. 
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Table 4. The observation log of the sample stars The right ascension (RA), decUnation (Dec) and date of the ISAAC observations are indicated. 
The date, UT and integration time T of the observation are given. We have assigned a quaUty label (Q) to the ISAAC spectra: G-good quality; 
L-low flux levels (<0.5 Jy); J-large flux-level jump (>50%) between the 3.3 and 3.5 micron spectrum. The seeing at A=0.5 /im, estimated by 
the VLT Astronomical Site Monitor (ASM), is given at the time of the 3.3 and 3.5 fjm ISAAC observations respectively. The final column lists 
the spectral calibrator stars. These standard stars have been observed during the same night as the science target, unless when flagged with 
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AA /AO /O AA /I 

09/02/2004 


00:36 


6 


J 


1 1 T / A O y( 

1.17/0.84 


y Ori 


o 
8 


HD z44o04 


05:31:57.25 


■ 1 1 . 1 1 /I 
+ 11:1 /:41.4 


AA /AO /OAAd 

09/02/2004 


AI .OA 

01:39 


6 


G 


A OO / 1 AT 

0.88/1.07 


y Ori* 


9 


RY Ori 


f\c o r\c\ n/' 

05:32:09.96 


AO /I A .1 O 

-02:49:46.8 


A .1 / 1 A/OAA,1 

04/10/2004 


06:34 


18 


G 


A CT / A y C\ 

0.57/0.60 


T Tl~\ A A C f\ A 

HD 44594 


10 


T Ori 


05:35:50.35 


ACO O . 1 C 1 

-05:28:35.1 


AC /I 1 /OAA1 

05/1 1/2003 


07:39 


6 


G 


1.72/ 1.66 


y Ori* 


1 1 

11 


V380 On 


05:36:25.43 


-06:42:57.7 


AC /I 1 /OAAO 

05/11/2003 


AO.AA 

08:00 


6 


G 


1 elA 1 'i £lf\ 

1.64/ 1.60 


y Ori* 


12 


V5o6 On 


05:36:59.25 


A/C. AA. 1 /I 

-06:09:16.4 


AC /II /OAAT 

05/1 1/2003 


08:27 


18 


G 


1 AO / I T 1 

1.98/ 1.71 


y Ori* 


13 


BF Ori 


05:37:13.26 


AZ O C AA Z 

-06:35:00.6 


AO / I A/OAA.1 

08/10/2004 


06:51 


18 


G 


1.46/ 1.64 


7 Ori 


14 


HD 37357 


05:37:47.08 


AZ /I O O A O 

-06:42:30.3 


AO / I A/OAA.1 

08/10/2004 


07:24 


14 


G 


1.56/ 1.89 


7 Ori 


15 


N3sk 81 


05:38:09.30 


-06:49:17.0 


08/10/2004 


09:06 


14 


G 


1.69/1.86 


y Ori 


16 


HD 3741 1 


AC O O . 1 y( CI 

05:38:14.51 


AC.O C 111 

-05:25:13.3 


AO /I A /O AA A 

08/10/2004 


08:38 


10 


G 


1.31/1.18 


y Ori 


1 T 
1 / 


Haro 13a 


ACQ o . 1 o on 

Uj: js: Is.zU 


— I) /:U2:2u.u 


TO /I A/TAA/1 
26/10/21)04 


ACQ A 





Cj 


A /fT / A A A 

0.42/ U. 44 


HD JOOJS 


18 


V599 On 


US:js:38.dU 


AT . 1 ^ . /I A 

-0/:lo:4o.O 


A A / 1 A /T AA/1 

09/10/2004 


08:36 


12 


G 


A ^ 1 / f\ n £^ 
O.ol / 0. lb 


e Cen* 


19 


V350 Ori 


AC /1A 11 "70 

05:40:11.78 


AA /I O 1 1 A 

-09:42:11.4 


A A / 1 A /O A A .1 

09/10/2004 


05:05 


18 


G 


1.13/ 1.48 


e Cen 


20 


HD 37806 


AC A 1 .AO OA 

05:41:02.29 


AO . /IT . A A T 

-02:43:00.7 


OO / I A/OAA/1 

28/10/2004 


A/:. A£ 

06:46 


6 


G 


A /I O / A /I C 

0.48 / 0.45 


HD 38858 


21 


HD 250550 


AZT.AI .CA AA 

06:01:59.99 


1 1 ZT.I A. CZT T 

+16:30:56.7 


1 O /AI /O AAC 

12/01/2005 


04:29 


6 


G 


1 TO / 1 TC 

1.78/1.75 


y Ori 


22 


VY Mon 


06:31:06.90 


1 A O AC O 

+ 10:26:05.3 


1 O /A 1 /O AAC 

12/01/2005 


04:50 


4 


G 


O AO / O AA 

2.08/2.00 


y Ori 


23 


LkHa 215 


A/C.n.4 1 OA 

06:32:41.80 


1 1 A. AA. T O 

+ 10:09:33.6 


O 1 /A 1 /OAAC 

21/01/2005 


01:55 


8 


G 


A n/Z. 1 A TA 

0.76/0.70 


f Ori 










O^ /AO /OAAC 

26/02/2005 


02:08 


8 




1 OA / I /I O 

1.20/1.42 


p Cen 


24 


HD 259431 


AZ" 1 '5 AC 1 A 

06:33:05.19 


■ 1 A. 1 A. OA A 

+10:19:20.0 


1 /A 1 /OAAC 

21/01/2005 


02:54 


6 


G 


A OO /A TA 

0.88/0.79 


/ILep 










OC /AO /OAAC 

25/02/2005 


AO . A'^ 

03:42 


6 




A iCT / A CT 

0.67/0.57 


e Cen 


25 


R Mon 


r\/' OA AA OA 

06:39:09.89 


. AO A A 1 A A 

+08:44:10.0 


1 A /A 1 /OAAC 

10/01/2005 


05:32 


6 


G 


A AZ / A TO 

0.96/0.72 


y Ori 


26 


V590 Mon 


A/" /I A /i /I CI 

06:40:44.57 


. AA A O AO 1 

+09:48:02.1 


1 A /A 1 /OAAC 

10/01/2005 


03:46 


28 


G 


A TO / I A /I 

0.73/ 1.94 


y Ori 


27 


GU CMa 


AT - AI ./(A CI 

07:01:49.51 


1 1 . 1 O .AO O 

-11:18:03.3 


AI /AO/OAA/( 

01/03/2004 


03:48 


6 


J 


A ^ZT / A O y( 

0.66/0.84 


e Cen 


28 


HD 53367 


AT.A/f.'^C CO 

07:04:25.53 


1 A.OT. 1 C T 

-10:27:15.7 


1 O /AO /OAAC 

12/02/2005 


05:45 


6 


G 


1 AA / A TZ" 

1.00/0.76 


P Cir 










1 A /I 1 /OA AO 

19/11/2003 


08:15 


6 




A CA / A C 1 

0.50/0.51 


y On 


29 


NX Pup 


r\'~j 1 A o o o 

07:19:28.26 


A A O C 11 O 

-44:35:11.3 


A A /AO /OA A ,1 

09/02/2004 


05:04 


6 


G 


A O 1 / A O C 

0.81/0.85 


y Ori 


30 


HD 76534 


AO CC AO T 1 

08:55:08.71 


A O T7 CA A 

-43:27:59.9 


1 1 /A 1 /OAAC 

11/01/2005 


04:54 


18 


G 


O AO / I O C 

2.03/ 1.85 


V Pup 


31 


HD 85567 


AA.CA.'^O CA 

09:50:28.54 


^A.CO.AO A 

-60:58:03.0 


OA/A1 /OAAC 

20/01/2005 


09:06 


6 


G 


f\ y A 1 f\ CA 

0.64/0.59 


n Sco 


32 


HD 95881 


1 1 . A 1 . /CO 


T 1 . A. /I O A 

-/1:J0:48.4 


OO /A 1 /OAAC 

23/U1/20Uj 


05:34 


6 


G 


0.81/0.76 


^ Cen 


33 


HD 97048 


1 1 AO AO OO 

11:08:03.32 


T7 OA 1 "7 C 

-77:39:17.5 


1 O /AT /O A A .1 

1 8/07/2004 


01:01 


6 


G 


A C C / A i'O 

0.55/0.63 


e Cen 


34 


HD 100546 


11 OO OC /I/I 

1 1:33:25.44 


TA 11 .1 1 O 

-70:11:41.2 


OT /AT /OA A .1 

27/07/2004 


23:49 


6 


G 


1 O /I / 1 .11 

1.24/ 1.41 


K Eri 


35 


HD 101412 


11 OA A A Ay 

11:39:44.46 


/'A 1 A OT "I 

-60:10:27.7 


O O /AT /OA A /i 

28/07/2004 


00:20 


8 


G 


1 O O / 1 OA 

1.33/ 1.29 


/( Eri 


36 


T Cha 


1 1 CT. 1 O 

11:57:13.53 


TA.0 1 O 1 C 

-79:21:31.5 


1 O /AT /O AA /I 

18/07/2004 


01:26 


10 


J 


A OO /ATI 

0.88/0.71 


e Cen* 


^7 

J 1 


Loo 1 D 


1 ^-1 Q-n^ QS 

1 J. iy.uj.7o 


— DZ. JH". lU. i 




UJ. J7 


o 


T 


n 78 /n on 


6 L^eii 


38 


HD 130437 


14:50:50.26 


-60:17:10.3 


10/05/2004 


06:03 


8 


G 


0.93/0.84 


e Cen 


39 


SS73 44 


15:03:23.80 


-63:22:59.0 


05/09/2003 


01:08 


8 


G 


0.83/0.95 


y^Cet* 


40 


HD 132947 


15:04:56.05 


-63:07:52.6 


05/09/2003 


01:35 


8 


L+J 


0.88/0.93 


Cet* 


41 


HD 141569 


15:49:57.75 


-03:55:16.4 


22/08/2004 


23:15 


8 


J 


0.72/0.73 


(T Sco 


42 


HD 142666 


15:56:40.02 


-22:01:40.0 


23/08/2004 


00:53 


6 


J 


0.96/0.96 


cr Sco 


43 


HD 142527 


15:56:41.89 


-42:19:23.3 


05/09/2003 


02:05 


6 


G 


1.32/0.91 


y^Cet* 


44 


HR5999 


16:08:34.29 


-39:06:18.3 


22/07/2004 


03:46 


6 


J 


1.61/1.82 


HD 172910 


45 


HD 150193 


16:40:17.92 


-23:53:45.2 


10/05/2004 


06:50 


6 


G 


0.83/0.86 


/<• Aql 


46 


AKSco 


16:54:44.85 


-36:53:18.6 


10/05/2004 


07:15 


8 


G 


0.85/0.89 


*: Aql 


47 


V921 Sco 


16:59:06.90 


-42:42:08.0 


22/07/2004 


03:01 


6 


G 


1.76/1.82 


HD 172910 


48 


KKOph 


17:10:08.07 


-27:15:18.2 


22/07/2004 


03:26 


8 


J 


1.39/1.36 


HD 172910 


49 


HD 163296 


17:56:21.29 


-21:57:21.9 


19/09/2003 


02:07 


6 


G 


1.54/1.74 


/fCet 


50 


MWC297 


18:27:39.60 


-03:49:52.0 


19/09/2003 


02:27 


4 


G 


1.67/1.64 


/(Cet 


51 


VVSer 


18:28:47.87 


+00:08:39.6 


19/09/2003 


02:46 


6 


G 


1.42/1.32 


/(Cet 
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Table 4. — (Continued) 



Acke & van den Ancker: Diamonds in Herbig Ae/Be systems 



# 








Date 


UT 


T 


V 


Q PPi 1-1 o" 








hh:mm:ss 


dd:mm:ss 


dd/mm/yy 


hh:mm 


m 




arcsec 




52 


MWC 300 


18:29:25.69 


-06:04:37.3 


19/09/2003 


03:02 


6 


G 


1.42/1.45 


KCst 


53 


AS 310 


18:33:21.17 


-04:58:06.7 


19/09/2003 


03:26 


14 


L 


1.57/1.28 


/cCet 


54 


S CrA 


19:01:08.60 


-36:57:20.0 


19/09/2003 


03:48 


8 


G 


1.49/1.39 


K Cet 


55 


TY CrA 


19:01:40.83 


-36:52:33.9 


19/09/2003 


04:08 


8 


G 


1.14/1.04 


K Cet 


56 


RCrA 


19:01:53.65 


-36:57:07.6 


19/09/2003 


04:27 


6 


G 


0.83/0.91 


/f Cet 


57 


TCrA 


19:01:58.80 


-36:57:49.0 


26/09/2003 


01:38 


8 


G 


0.99/1.24 


BY Cap 


58 


HD 179218 


19:11:11.25 


+15:47:15.6 


01/10/2003 


02:01 


6 


G 


0.63/0.66 


HD9562 


59 


HBC 687 


19:29:00.80 


+09:38:43.0 


26/09/2003 


02:08 


20 


L 


1.11/1.12 


HD 9562 


60 


V1295 Aql 


20:03:02.52 


+05:44:16.7 


19/09/2003 


04:46 


6 


G 


0.78/0.90 


/cCet 
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Table 5. The dete cted features in the 3 micron region of the sample stars. The sample stars have been classified in Meeus groups I and II 
<Meeus et alJl200lh . Targets with Ay exceeding 4 are classified as embedded sources (E). Sources with unknown or uncertain classification 
are flagged with "?". The measured line fluxes and peak-over-continuum ratios of four spectral features are given: the Pfund 6 emission line 
at 3.296 yum, the PAH 3.3 ^m feature and the nanodiamond features at 3.43 and 3.53 ^im. Tentative detections (line flux between l-3cr) are 
indicated with a. The flag b indicates detected features of unclear origin, either linked to data reduction artefacts, nanodiamonds or an unknown 
carrier. 



# 


Target 


Group 


Pf 6 


PAH 


NANO 


NANO 


Pf S 


PAH 


NANO 


X T A X T/^ 

NANO 








3.296 yum 


3.3 yum 


3.43 ^im 


3.53 yum 


3.296 yum 


3.3 yum 


3.43 yum 


3.53 yum 










fine flux/ 10" 


'° Wm ^ 






peak-over-continuum ratio 




1 


UX Ori 


II 


0.7 ± 0.2 


< 14 


< 16 


< 2.9 


1.18 ± 0.05 


< 1.05 


< 


1.05 


< 1.01 


2 


HD 34282 


I 


< 1.1 


24 ± 7 


< 17 


< 4.7 


< 1.06 


1.22 ± 0.07 


< 


1.06 


< 1.02 


3 


V346 Ori 


I 


< 1.4 


17 ± 8" 


< 23 


< 8.4 


< 1.10 


1.2 ± 0.1 


< 


1.10 


< 1.05 


4 


CO Ori 


II 


< 0.4 


< 6.9 


< 8.2 


1.2 ± 0.7"-* 


< 1.05 


< 1.05 


< 


1.15 


1.02 ± 0.02 


5 


HD 35929 


II 


0.4 ± 0.4" 


< 12 


< 14 


< 2.7 


1.04 ± 0.04 


< 1.04 


< 


1.04 


< 1.01 


6 


StHa 41 


II? 


< 0.32 


< 5.8 


< 10 


< 4.2 


< 1.05 


< 1.05 


< 


1.15 


< 1.05 


7 


HK Ori 


II 


< 1.1 


< 31 


< 35 


< 9.0 


< 1.08 


< 1.08 


< 


1.08 


< 1.03 


8 


HD 244604 


II 


< 1.5 


< 27 


< 31 


< 8.0 


< 1.10 


< 1.10 


< 


1.10 


< 1.04 


9 


RY Ori 


II 


< 0.25 


< 4.8 


< 5.7 


< 1.1 


< 1.06 


< 1.06 


< 


1.06 


< 1.02 


10 


T Ori 


II? 


< 3.7 


< 62 


< 69 


< 12 


< 1.06 


< 1.06 


< 


1.06 


< 1.01 


11 


V380 Ori 


II? 


5 + 1 


< 86 


< 110 


< 22 


1.27 ± 0.08 


< 1.07 


< 


1.08 


< 1.02 


12 


V586 Ori 


II 


0.25 ± 0.09° 


< 8.7 


< 10 


< 3.0 


1.13 ± 0.05 


< 1.05 


< 


1.05 


< 1.02 


13 


BP Ori 


II 


0.6 ± 0.2 


< 7.6 


< 8.8 


< 3.5 


1.15 ± 0.05 


< 1.05 


< 


1.05 


< 1.03 


14 


HD 37357 


II 


2.1 ± 0.4 


< 6.2 


< 7.1 


< 3.9 


1.15 ± 0.03 


< 1.03 


< 


1.03 


< 1.02 


15 


N3sk 81 


E 


< 0.58 


< 9.7 


< 11 


< 3.1 


< 1.04 


< 1.04 


< 


1.04 


< 1.02 


16 


HD 37411 


II 


< 0.65 


< 10 


< 12 


< 4.4 


< 1.04 


< 1.04 


< 


1.04 


< 1.02 


17 


Haro 13a 


I 


< 1.9 


< 34 


< 44 


< 7.3 


< 1.03 


< 1.03 


< 


1.03 


< 1.01 


18 


V599 Ori 


E 


< 0.24 


5 ± 1 


< 4.2 


< 0.95 


< 1.04 


1.15 ± 0.04 


< 


1.04 


< 1.01 


19 


V350 Ori 


II 


< 0.063 


< 1.1 


< 1.2 


< 0.37 


< 1.06 


< 1.06 


< 


1.06 


< 1.02 


20 


HD 37806 


II 


1.3 ± 0.3 


< 17 


< 20 


< 2.3 


1.14 ± 0.03 


< 1.03 


< 


1.03 


< 1.01 


21 


HD 250550 


I 


< 3.3 


< 59 


< 69 


< 8.3 


< 1.16 


< 1.16 


< 


1.16 


< 1.03 


22 


VY Mon 


E 


< 31 


< 530 


< 670 


< 38 


< 1.17 


< 1.17 


< 


1.17 


< 1.01 


23 


LkHtt 215 


I 


< 0.40 


< 7.3 


< 8.8 


< 1.7 


< 1.07 


< 1.07 


< 


1.07 


< 1.02 


24 


HD 259431 


I 


0.6 ± 0.2" 


< 32 


< 37 


< 4.5 


1.16 ± 0.06 


< 1.06 


< 


1.06 


< 1.01 


25 


RMon 


I 


< 13 


< 260 


< 350 


28 ± lO"* 


< 1.08 


< 1.08 


< 


1.08 


1.03 ± 0.01 


26 


V590 Mon 


7 


< 1.9 


7 ± 6" 


< 34 


2.4 ± 0.8"-* 


< 1.16 


1.2 ± 0.2 


< 


1.16 


1.05 ± 0.02 


27 


GU CMa 


II 


4.9 ± 0.8 


< 43 


< 48 


< 12 


1.46 ± 0.09 


< 1.08 


< 


1.08 


< 1.03 


28 


HD 53367 


7 


< 1.3 


< 20 


< 24 


< 5.1 


< 1.11 


< 1.11 


< 


1.11 


< 1.03 


29 


NX Pup 


II 


< 3.4 


< 67 


< 75 


< 15 


< 1.04 


< 1.04 


< 


1.04 


< 1.01 


30 


HD 76534 


7 


< 0.27 


< 4.5 


< 5.2 


< 0.59 


< 1.19 


< 1.19 


< 


1.19 


< 1.03 


31 


HD 85567 


II 


< 1.2 


< 19 


< 21 


< 3.0 


< 1.06 


< 1.06 


< 


1.06 


< 1.01 


32 


HD 95881 


II 


< 3.1 


< 49 


< 56 


< 9.9 


< 1.10 


< 1.10 


< 


1.10 


< 1.02 


33 


HD 97048 


I 


< 1.5 


65 ± 20 


200 ± 20 


486 ± 5 


< 1.10 


1.4 ± 0.1 


2.0 


± 0.1 


4.3 ± 0.2 


34 


HD 100546 


I 


9 ± 2 


240 ± 40 


< 94 


13 ± 10"* 


1.33 ± 0.06 


1.34 ± 0.06 


< 


1.06 


1.05 ± 0.04 


35 


HD 101412 


II 


< 1.3 


13 ± 10" 


< 25 


< 7.5 


< 1.06 


1 AO 1 n f\£^ 

1.08 ± 0.06 


< 


1.06 


< 1.03 


36 


T Cha 


II 


< 0.85 


< 16 


< 17 


< 4.0 


< 1.07 


< 1.07 


< 


1.07 


< 1.02 


37 


LSS 3027 B 


7 


< 0.46 


< 8.7 


< 8.8 


< 3.6 


< 1.73 


< 1.73 


< 


1.73 


< 1.39 


DO 


rlU 1 jU4j / 




l.o ± u.z 


< D.J 


< /.i 


< J.O 


1.04 ± U.Uo 


< 1.06 


< 


1.06 


^ 1 nA 
< i.uo 


39 


SS73 44 


I 


<0.89 


< 16 


< 17 


< 8.8 


< 1.05 


< 1.05 


< 


1.05 


< 1.04 


40 


HD 132947 


II 


< 0.68 


< 8.6 


<9.1 


< 8.1 


< 1.14 


< 1.14 


< 


1.14 


< 1.18 


41 


HD 141569 


II 


<0.87 


<37 


<38 


<5.9 


< 1.20 


< 1.20 


< 


1.20 


< 1.05 


42 


HD 142666 


II 


<0.72 


< 20 


< 21 


< 5.8 


< 1.04 


< 1.04 


< 


1.04 


< 1.02 


43 


HD 142527 


I 


32 ±5 


< 110 


< 120 


44 + 9* 


1.27 ±0.04 


< 1.04 


< 


1.04 


1.05 ±0.01 


44 


HR 5999 


II 


5±3" 


< 110 


< 150 


<55 


1.08 ±0.04 


< 1.04 


< 


1.04 


< 1.02 


45 


HD 150193 


II 


3.4 ±0.3 


< 13 


< 15 


<6.3 


1.15 ±0.01 


< 1.01 


< 


1.01 


< 1.01 


46 


AK Sco 


II 


<0.79 


< 13 


< 15 


2.5 ± 0.9"-* 


< 1.03 


< 1.03 


< 


1.03 


1.03 ±0.01 


47 


V921 Sco 


E 


86 ±5 


< 160 


< 190 


103 ± 30* 


1.42 ±0.03 


< 1.03 


< 


1.08 


1.04 ±0.01 


48 


KK Oph 


II 


<2.8 


<47 


<66 


< 19 


< 1.05 


< 1.05 


< 


1.05 


<1.02 


49 


HD 163296 


II 


18 ± 3 


< 120 


< 130 


23 ±7* 


1.16 ±0.03 


< 1.03 


< 


1.03 


1.03 ±0.01 


50 


MWC 297 


E 


240 + 20 


<910 


< 1100 


535 ± 100 


1.63 ±0.05 


< 1.04 


< 


1.04 


1.06 ± 0.01 


51 


VV Ser 


II 


2.1 ±0.4 


<23 


<25 


<9.0 


1.13 ±0.03 


< 1.03 


< 


1.03 


<1.01 


52 


MWC 300 


E 


9± 1 


< 82 


<97 


50 ±9* 


1.13 ±0.02 


< 1.02 


< 


1.02 


1.03 ±0.01 


53 


AS 310 


I 


<0.83 


11 ± 3 


<9.5 


<4.0 


< 1.34 


2.1 ±0.5 


< 


1.34 


< 1.22 
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# 


Target 


Group 


Pf 5 


PAH 


NANO 


NANO 


Pf 5 


PAH 


NANO 


NANO 








3.296 


3.3 yum 


3.43 yum 


3.53 yum 


3.296 yum 


3.3 yum 


3.43 yum 


3.53 yum 










line flux/ 10" 


Wm-2 






peak-over-continuum ratio 




54 


SCrA 


7 


3.9 ± 0.4 


<32 


<35 


10 + 3" 


1.20 ±0.03 


< 1.03 


< 1.03 


1.03 ±0.01 


55 


TY CrA 


7 


<0.75 


40 ±4 


< 12 


<4.7 


< 1.03 


1.28 ± 0.03 


< 1.03 


< 1.02 


56 


RCrA 


1 


<26 


<470 


<550 


290 + 50* 


< 1.02 


< 1.02 


< 1.02 


1.04 ±0.01 


57 


TCrA 


II? 


<2.5 


<44 


<96 


7 ± 1* 


< 1.06 


< 1.06 


< 1.19 


1.04 ±0.01 


58 


HD 179218 


I 


<3.1 


160 ± 20 


<47 


<6 


< 1.03 


1.28 ± 0.03 


< 1.03 


< 1.04 


59 


HBC 687 


I? 


<0.14 


1.7± 0.7" 


<2.8 


< 1.0 


< 1.07 


1.14 ±0.07 


< 1.07 


< 1.03 


60 


VI 295 Aql 


II 


4.8 ± 0.6 


<40 


<44 


< 16 


1.18 ±0.02 


< 1.02 


< 1.02 


< 1.01 
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